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ABSTRAGT 

We study dynamical capture binary black hole-neutron star (BH-NS) mergers focusing on the 
effects of the neutron star spin. These events may arise in dense stellar regions, such as globular 
clusters, where the majority of neutron stars are expected to be rapidly rotating. We initialize the 
BH-NS systems with positions and velocities corresponding to marginally unbound Newtonian orbits, 
and evolve them using general-relativistic hydro dynamical simulations. We find that even moderate 
spins can significantly increase the amount of mass in unbound material. In some of the more extreme 
cases, there can be up to a third of a solar mass in unbound matter. Similarly, large amounts of 
tidally stripped material can remain bound and eventually accrete onto the BH—as much as a tenth 
of a solar mass in some cases. These simulations demonstrate that it is important to treat neutron 
star spin in order to make reliable predictions of the gravitational wave and electromagnetic transient 
signals accompanying these sources. 

Subject headings: black hole physics—gamma-ray burst: general—gravitation—gravitational waves— 
stars: neutron 


1. INTRODUCTION 

Next generation ground-based g ravitational wave 
(GW ) detectors such as aLIGO (|Abramovici et al.l 
are expected to reach design sensitivity within 
the next few years. Among their most promis¬ 
ing sources are mergers of compact objects (COs), 
including black hole-neutron star (BH-NS) bina¬ 
ries. BH-NS mergers are also propos ed short- 
hard gamma-ray burst (sGRB) engines (e.g. iMeszaro^ 
(1200^ ). and may power ot her electromagnetic (EM) 
transients either preceding (iHansen fc LvutikovI l2QQll: 
iMcWilliarns fc LevinI l2QllI iPaschalidis et al.l l2Q13f ) or 
following ( Metzger fc BergCT 120121) the merger. These 
EM counterparts to GWs could be observed by current 
and future wi de-held telescopes such as PTE (|Ruull2QQ9l) 
and LSST (|LSST Dark Energy Science CollaborationI 

M2). 

Extracting maximum information from such “mul¬ 
timessenger” observations requires careful modeling of 
BH-NS mergers. Several studies of quasicircular BH-NS 
inspirals using nume rical relativity s imula t ions have been 
perfor m ed, see, e.g..lChawla et al. (I2Q1QD: [Etienne et 4] 
(1 201 2D : Ik vutoku et al .1 (|2Q13D : Tanaka et al.l (j2Q14D : 
lEoucart et al.l (|2Q14D . While quasicircular binaries may 
dominate the global rates of BH-NS encounters in 
the Universe, recent calculations (iKocsis fc LevinI I2Q12I: 
iLee et all I2Q1QI: iSamsing et al.l I2Q14D suggest that in 
dense stellar regions, such as galactic nuclei and globular 
clusters (GGs), CO binaries can form through dynami¬ 
cal capture and merge with non-negligible eccentricities. 
Compared to quasi-circular inspirals, these emit more 
GW energy in the high luminosity, strong-held regime of 
general relativity, and small changes in the energy of the 
binary at each pericenter passage can lead to relatively 
large changes in the time between GW bursts—the lead¬ 
ing order GW observable. Hence these systems could 
be excellent laboratories to test gravity and measure the 


internal structure of a NS (insofar as this ahects the en¬ 
ergy of the orbit, e.g. tidal excitation of f-modes). Rates 
of these events are highly uncertain, but have been es¬ 
timated to be up to ^ 100 yr“^ Gpc“^. To realize this 
rich potential to learn about the Universe from eccentric 
mergers, it is irrelevant what their rates are compared 
to quasi-circular inspiral, only that eccentric mergers oc¬ 
cur frequently enough that some events could be observ¬ 
able by aLIGO within its lifetime. However, new detec¬ 
tion pipelines would be needed that are better adapted 
to the repeated-burst nature of eccent ric GW mergep 
for aLIGO to efficiently detect them (jTai et al.l I2014D . 
Eor more discussion of rates, distinguishing fe atures, and 
detect ion issues for eccentric encounters, see lEast et ^ 
(HoH) and references therein. 

Motivate d by the above, [Stephens et al.l 
lEast et, aH (l2fll2aD performed fully general-relativistic 
hydrodynamical (GR-HD) simulations of dynamical- 
capture BH-NS mergers. These studies explored the ef¬ 
fects of impact parameter, BH spin, and NS equation 
of state (EOS) on GW emission and post-merger BH 
disk and ejecta masses. Here, we expand upon this work 
by including the effects of NS spin. To date the only 
simulations including spin ning N S s focu s ed on quasicir¬ 
cular NS-NS mergers, e.g. iTichvI (j2QllD : iBernuzzi et al.l 
(IMl . demonstrating that even moderate spins can af¬ 
fect the dynamics. 

NS spin has two main effects: (1) it modifies the star’s 
structure, making it less gravitationally bound; (2) it 
changes the orbital dynamics, e.g., by shifting the effec¬ 
tive innermost stable orbit (ISO). This can impact not 
only the GWs from CO mergers, but also the amount 
of matter forming the BH accretion disk that putatively 
powers a sGRB, and the amount of unbound matter that 
powers other EM transients, such as kilonovae. There 
may also be effects on pre-merger EM signals since the 
NS spin determines the light-cylinder radius, and hence 




























































2 


East, Paschalidis & Pretorius 


the orbital separation at which unipolar induction turns 


Spin effects on the NS structure cannot be neglected if 
the NS spin period P is O(nis). Furthermore, for com¬ 
parable mass BH-NSs near the tidal disruption radius, 
NS spin effects on the orbit will be non-negligible when 
P is similar to the BH-NS encounter timescale di ic 
1201 111 . For example, a BH-NS eccentric encounter with 
mass ratio q = Mbu/M^s = 4 (as studied here) near 
a periapse of = lOM has an interaction timescale 
of tint — ^ 1.0(MNs/l-4M©)ms (M is the 

system’s total mass, and we use geometric units with 
G = c = 1 throughout). 

NSs in field BH-NS binaries may not commonly have 
P = 0(ms) near merger. However, there are two reasons 
to think that the opposite may hold for dynamical cap¬ 
ture BH-NS mergers occurring in GCs: the pulsar spin 
period distribution in Galactic GCs peaks in the mil¬ 
liseconds, and millisecond pulsars (MSPs) have longer 
inferred magnetic dipole spin-down timescales. 

Of the 144 currently known pulsars in Galactic GCs, 
^ 83% have periods less than 10 ms, ^ 55% less than 5 
ms, and ~ 12% have periods less than 2.5 rnsQ. This set 
includes PSR-J1748-2446ad—the fastest-spinning pn l- 
sar known, with P = 1.396 ms (jHessels et al.l 1200^ . 
The theoretical explanation for this skew toward short 
periods is that GCs fav or the formation of l ow-mass 
X-ray binaries (LMXB) (jVerbnnt fc HutI 119871 ). which 
are thought to spin up the NS to m s periods through 
mass and angular momentnni trans fer (jAlpar et al.lll982l : 
iRadhakrishnan Sz Srinivasanlll982f ). 

Assuming that pulsar spin-down is predominantly due 
to magnetic dipole emission, both the magnetic field 
strength (B) and spin-down timescale (tsd) can be com¬ 
puted from observation s of P and its time derivative P. 
For B the relation is (jBhattacharva fc van den Henvell 

[Mil) 


B 


1.6 X 10® G 


2.5 ms 10-20 


1/2 


which for known GC MSPs with P > 0 gives typi¬ 
cal values B ^ 10^—lO^G. Observations of X-ray os¬ 
cillations of accretion-powered MSPs in LMXBs, and 
pulsar recycli ng theory imply B in the range 3 x 
10 ^— 3 X IQ^G (iLamb fc Yul2QQ5[) . For Gd the expression 
is (|Zhang fc Meszarosll2QQl ) 


tsd ^ 4 Gyr 


1045g 


B 


2.5 ms 


A 

10 km / 


3 X 10^ G 

-6 


-2 


( 1 ) 


where the NS moment of inertia is 1. Even neglecting the 
possibility of magnetic field decay, e.g. a pulsar with a 
magnetic field of 3 x lO^G and initial P = 2.5ms will take 
roughly a Hubble time for P to double. Give n the long 
spin-d own timescale of MSPs and the results of iLee et al.l 
(l20T0h . which suggest that in GCs there could be 40 BH- 
NS collisions per Gyr per Milky Way-equivalent galaxy. 


www.naic.edu/~pfreire/GCpsr.html 


it is at least conceivable that some of these eccentric BH- 
NS collisions take place with millisecond NSs. 

With this motivation, here we focus on eccentric BH- 
NS mergers (with initial conditions corresponding to a 
marginally unbound Newtonian orbit) and explore NS 
spin effects. We show that even moderate spins can 
significantly impact the outcome, both in terms of the 
GWs, and amounts of tidally stripped bound and un¬ 
bound matter. The remainder of the paper is as follows: 
in Sec.[2]we describe our initial data and numerical meth¬ 
ods. In Sec. [3] we present our simulation results and 
discuss the impact of NS spin on gravitational and EM 
signatures. We summarize in Sec. |4] and discuss future 
work. 

2. NUMERICAL APPROACH 

We perform GR-HD simulati ons of BHs mergin g with 
rotating NSs using the code of lEast et al.l (j2Q12bf ). The 
field equations are solved in the generalized-harmonic 
formulation, using finite differences, while the hydrody¬ 
namics are evolved using t he same hig h -resolu tion shock¬ 
capturing techniques as in lEast et all (j2Q12aI ) . 

To construct initial data, we solve the c onstraint equa¬ 
tions using the code of lEast et al.l (j2Q12c[ ). specifying the 
free-data as a superposition of a non-spinning BH with an 
equilibrium, uni formly rotating NS , which we generate 
using the code of iCook et al.l (jl994allbl ). For t he NS EOS, 
we ad opt the HB piece-wise polytrope from iRead et al l 
(j2QQ9f ). and include a thermal component Pth = 0.5ethP 
allowing for shock heating. 

Fixing the NS gravitational mass to 1.35M0, we con¬ 
sider NSs with dimensionless spins uns = J/^^s = 0? 
0.1, 0.2, 0.3, 0.4, 0.5, and 0.756, having correspond¬ 
ing compactions (mass-to-equatorial-radius) C = 0.172, 
0.171, 0.169, 0.166, 0.161, 0.154, and 0.12. The ranges 
of T/|1T| (kinetic-to-gravitational-potential-energy ra¬ 
tio) and P in our spinning NS models are [0.003,0.12] 
and [5.25,1.00] ms. The fastest spinning NS considered 
has a polar-to-equatorial-radius ratio of r^o/^eg = 0.55, 
near the mass-shedding limit of r^o/^eg = 0.543. We also 
vary Vp/M G [5, 8]. 

We consider systems with q = A. The two COs are ini¬ 
tially placed at a separation of d = 50M (^ 500km), with 
positions and velocities corresponding to a marginally 
unbound Newtonian orbit labeled by Vp. For this initial 
study, we only consider cases where the spin is aligned 
or anti-aligned with the orbital angular momentum (the 
latter indicated by uns < 0)- 

The simulations utilize seven levels of adaptive mesh 
refinement that are dynamically adjusted based on the 
estimated truncation error of the metric. Most simula¬ 
tions are performed using a base-level resolution with 
193^ points, and finest-level resolution with approxi¬ 
mately 75 (130) points covering the (non-spinning) NS 
(BH) diameter, respectively. For Vp/M = 6, uns = 0.756 
we also perform simulations at 2/3 and 4/3 x the reso¬ 
lution, to establish convergence and estimate truncation 
error. In Fig. [T] we demonstrate the convergence in the 
GW emission. 

3. RESULTS AND DISCUSSION 
3.1. Simple estimates 

During eccentric encounters between NSs and BHs, for 
NS tidal disruption to form a substantial disk and un- 
















































Dynamical Capture Mergers of BHs with spinning NSs 


3 




Fig. 1.— Convergence of the GW emission for the Vp/M = 6 , 
dNS = 0.756 case. The top panel shows the real part of the i = m = 
2 mode of the Newman-Penrose scalar T 4 at three resolutions, and 
the bottom panel the differences in this quantity with resolution, 
scaled assuming second-order convergence. 

bind a non-negligible amount of matter, it must occur 
outside the ISO radius (nso = 4 Mbh foi* a marginally 
unbound test particle about a non-spinning BH). In addi¬ 
tion to shifting the effective ISO, spin makes the NS less 
self-bound, and thus alters the tidal disruption radius. 
Equating the sum of the tidal and centrifugal accelera¬ 
tions to the gravitational acceleration on the NS surface 
yields 


le-2 



where we replaced the NS angular frequency with Cl = 
Jjl — uns^ns/^’ ^ ~ 2/Mns^ns/^’ with / an 

order-unity constant that depends on the NS structure. 
Here, a^s = (4/2/250)1/2 ~ 0.8(//0.7)(C/0.12)-i/2 is 
the mass-shedding limit spin parameter. Equation m 
shows that the closer uns is to ams, the larger the tidal 
disruption radius. It also suggests that for sufficiently 
fast rotators, the tidal disruption radius can be outside 
the ISO even for large something which is not true 
for non-spinning NSs unless the BH has near-extremal 
spin. Additionally, as the prograde NS spin increases, 
the effective ISO decreases. Therefore, we expect more 
massive disks and more unbound material following tidal 
disruption outside the ISO with increasing uns- 

3.2. Dynamics and Gravitational Waves 

Eor the cases considered here, those with Vp < 6.5 
merge on the initial encounter, while those with Vp > 7.5 
go back out on an elliptic orbit after fly-by. Eigure [2] 
plots the dominant contribution to the GW signal for 
Vp/M = 6, 7, and 8; see also Table [TJ Near the critical 
Vp — below (above) where a merger (fly-by) occurs — 
there are large differences in the dynamics that have a 
noticeable impact on the GW signal and tidally stripped 
matter. This is evident here with the VpjM = 7 case, 
where there is either partial tidal disruption followed by a 
merger on a second encounter (for uns = 0.4 and 0.5), or 
complete tidal disruption/merger on the first encounter 
(for the other spins). This is illustrated in Eig. [3] for 



Fig. 2. — Plots of the real part of the I = m = 2 mode of the 
scalar T 4 . The top and middle panels show GWs from simulations 
with Vp/M = 6 and 7, respectively. The bottom panel plots the 
GWs for Vp/M = 8 ffy-by cases. 


= 0.2, 0.5, as well as the bottom-right panel of Eig.H] 
where it can be seen that the NS for uns = 0.4 (0.5) loses 
^ 10% (^ 20%) of its mass in the initial encounter. 

Though it is difficult to disentangle the nonlinear dy¬ 
namics occurring here, we suggest the following expla¬ 
nation for this non-monotonic behavior as a function of 
uns for TpjM — 7. The uns = 0.756 NS is the least 
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bound and for this case complete tidal disruption occurs 
on the initial encounter. As the NS is tidally stretched, 
GW emission effectively shuts off (middle panel Fig. [2]), 
and matter begins to accrete onto the BH (bottom panel 
of Fig. HI). For the next two lower spins, the material 
is more tightly bound, and only partial disruption oc¬ 
curs, with some of the material promptly accreting onto 
the BH. However, as the NS spin decreases, the effective 
ISO radius increases, and for uns ^0-3 the core of the 
NS crosses the ISO, resulting in immediate merger. The 
= 0.4 and 0.5 cases lose enough orbital energy and 
angular momentum during the first encounter that they 
merge on the second. 

Note that the simple estimate of Equation (|2]) implies 
that since at Vp/M = 7 the uns = 0.756 case is com¬ 
pletely disrupted, while some of the lower spin cases are 
partially disrupted, Vp/M = 8,aNS = 0.756 should also 
be disrupted, given how close this spin is to the break up 
value. That this does not happen shows that this crude 
estimate significantly underestimates the self-binding of 
high spin stars. 


3.3. Post-merger Matter Distribution 

In Table [T] we list the amount of bound and unbound 
mass exterior to the BH shortly following merger. For 
rpjM — 5,6 the bound mass is only a weak function of 
uns, with the notable exception of uns = 0.756, r^/M = 
6. By contrast, near the critical impact parameter 
{rpjM = 7) there is over an order of magnitude vari¬ 
ation in bound material as a function of NS spin. 

The amount of bound rest-mass that forms a BH accre¬ 
tion disk is 0.01-0.15M0 in our set. If these disks power 
sGRBs on timescales of ^ 0.2s, the accretion rates will 
be ^ 0.05-0.75 Mq s“^, i.e., c onsistent with magneto - 
hydrodynamic BH-NS studies (jPaschalidis et al.ll2014f ). 
Assuming a 1% conversion efficiency of accretion power 
to jet luminosity, these accretion rates imply luminosi¬ 
ties of lO^^-lO^^erg s“^—consistent with characteristic 
sGRB luminosities. 

The top two panels in Fig. |4] show plots of the asymp¬ 
totic velocity distribution of the un boun d matter for 
rpjM = 6, 7. As anticipated (see Sec. (j3.1j) ). the general 
trend is that increasing uns increases both the amount 
and average asymptotic velocity of unbound material. 
This is also seen in Table [T] where we list these quanti¬ 
ties for all cases. The bottom-left panel of Fig. |4] also 
demonstrates that including spin increases the amount 
of unbound material by an order of magnitude or more 
for the cases considered. 

The increase in total rest-mass Mo,u and velocity v of 
the unbound material with increasing uns can strongly 
impact potential kilonovae signatures from such mergers. 
These arise when neutron-rich ejecta produce heavy el¬ 
ements through t he r-process than then undergo fission , 
emitting photons jET h Paczvnskilll99§: lKnlkarnill2QQ5[ ). 
R.ecentlv iBarnes fc KasenI (j2Q13[ ) have shown that the 
opacities in r-process ejecta will likely be dominated by 
lanthanides, giving rise times of 


Vak « 0.25(Mo,VlO“" Mo)^/2(^/0.3c)-i/ 2 d 
with peak luminosities of 


L fa 2 X 10"^^ 


Mo.u 
10-2 M© 



1/2 


erg s 


(3) 


for typical values found here. In some cases, opacities 
an order of magnitude lower tha n those used above may 
be justified (jMetzger et al.l[2QT^ . Using Eq. m we esti¬ 
mate the luminosity from potential kilonovae in Table [TJ 
Apart from the fact that NS spin can make the difference 
in whether there will be a kilonova at all for rpjM = 5, 
for rpjM = 6 and rpjM = 7 spin affects L by an order 
of magnitude in our set. For rpjM = 6 even a moder¬ 
ate uns = 0.1 increases L by a factor of 4 compared to 
= 0. iBarnes h KasenI (j2Q13f ) predict that a kilonova 
luminosity of ^ lO^^erg s“^ corresponds to an r-band 
magnitude of 23.5 mag at 200 Mpc (near the edge of the 
aLIGO volume), above the planned LSST survey sensi¬ 
tivity of 24.5 mag. Thus, differences in luminosity by 
factors of a few could be discernible. 

Ejecta will also sweep the interstellar medium (ISM) 
producing radio waves. These will peak on t imescales of 
weeks with brightness (jNakar fc Piranl[2011f ) 

F{uohs)^0.6{E]^in/l0^^ erg)(no/0.1 cm“^)'^/^ (4) 

(n/0.3c)^^/^(z^obs/GHz)“^/^(d/100 Mpc)“^ mJy 

for an observation frequency z^obs at a distance d, 
and using no ^ 0 . 1 cm ~^ as the density for GG 
cores (|R.osswog et al.|[20l3 ). Estimating the kinetic en¬ 
ergy and the mass-averaged velocity in the ejecta, we 
show F(z/obs) via Eq. (|4]) in Table [TJ For rpjM = 6 
{rpjM — 7) F(z/obs) varies by 3 (2) orders of magnitude 
over our set of spins. 

Finally, it has been suggested that ejecta from merg¬ 
ers involving NSs may make a non-negligible contri- 
bntion to the overall abundance of r-process elem ents 
(|Lattimer fc Schramm! 11974 lR.osswog et al.l Il998f ) . In 
particular, dynamical-capture binaries, which can form 
and merge on shorter timescales, may be favored 
over field binaries in explai ning abundances in cpbon - 
enhanced metal-poor stars (|Rnmirez-Rniz et ahl l2014f ) . 
The average galactic production o f these ele ments is 
estimated to be ^ yr“ ^ (|Qianl l2QQQh . Mak¬ 

ing the limiting assumption that all r-process mate¬ 
rial comes from extreme BH-NS mergers cases like the 
rp = 7, ttNS = 0-T6 with Mo,u ~ 0.3Mq caps these 
extreme events at 3 x 10“^ per galaxy (sim¬ 

ilar to predicted rate s for primordial BH-NS merg¬ 
ers (jAbadie et al.ll2QTQh ). 

4. CONCLUSIONS 

We have demonstrated using GR-HD simulations of 
dynamical capture BH-NS mergers that even moderate 
values of NS spin can significantly increase the mean ve¬ 
locity and amount of unbound material (to as much as 
O.3M0 for extreme spins). This could lead to signifi¬ 
cantly brighter transients, including kilonovae a factor 
of a few brighter, and radio wave emission from inter¬ 
action with the ISM an order of magnitude or more 
brighter. For comparison, simulations of quasicircular 
BH-NS mergers with nonspinning NSs typically find 
ejecta velocities ^ 0.2-0.3c, comparable, though some¬ 
what smaller than found here, but only find similar 
amounts of ejected material f or cases with smaller mass- 
ratios and/or high BH spin (jKvntokn et al.ll2015f ). We 
also find that the NS spin can alter the amount of bound 
matter that, following tidal disruption, remains to form 
an accretion disk that may power a sGRB. Depending 
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Fig. 3. — Equatorial density snapshots. Top row (r^/M = 7, uns = 0 - 5 ) from left to right: the NS survives the first encounter (first 
and second panels), it is completely tidally disrupted during the second encounter (third panel), the bulk of the matter outside the BH 
is unbound (fourth panel). Bottom row (r^/M = 7, uns = 0-2) from left to right: the NS is tidally disrupted during the first encounter 
(first panel), a tidal tail forms ejecting some matter to infinity (second panel), an accretion disk develops outside the BH (third and fourth 
panels). The scale can be inferred from the size of the BH (^bh ~ 16 km). 






Fig. 4. — Top: distribution of the asymptotic velocity of unbound rest-mass, binned in increments of 0.05c, and computed 10 ms 
post-merger for r^/M = 6 (left) and r^/M = 7 (right) and various spins. Bottom: total unbound mass as a function of NS spin (left) and 
rest-mass outside the BH versus time for Vp/M = 7 and various spins (right). Note that for the lowest point in the bottom-left panel, no 
unbound matter was found in the simulation (indicated by an arrow). 
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5.0 
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0.81 
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0.94 
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0.19 

0.3 

1.3 

0.05 

0.51 

5.0 

0.75 

0.55 

0.81 

4.97 

1.12 

1.51 

0.42 

3.0 

2.9 

4.5 

0.51 

6.0 

-0.40 

0.55 

0.92 

5.66 

1.01 

0.50 

0.22 

0.3 

1.2 

0.08 

0.51 

6.0 

0.00 

0.56 

1.13 

6.86 

1.15 

0.02 

0.18 

0.007 

0.2 

0.001 

0.52 

6.0 

0.10 

0.57 

1.16 

7.47 

1.07 

0.30 

0.19 

0.1 

0.9 

0.02 

0.53 

6.0 

0.20 

0.57 

1.21 

7.30 

1.10 

0.21 

0.20 

0.1 

0.7 

0.02 

0.54 

6.0 

0.40 

0.58 

1.31 

8.03 

0.91 

1.99 

0.32 

2.3 

2.9 

1.6 

0.53 

6.0 

0.75 

0.60 

1.04(1.23)^ 

6.85(7.92) 

2.40(2.23) 

14.39(14.01) 

0.346(0.350) 

19.9 

8.2 

17.9 

0.491 (0.492) 

6.5 

0.40 

0.60 

1.51 

9.91 

4.39 

9.08 

0.28 

8.0 

5.8 

3.8 

0.50 

7.0 

-0.40 

0.63 

1.51 

9.38 

1.17 

0.35 

0.23 

0.2 

1.0 

0.06 

0.53 

7.0 

0.00 

0.61 

1.72 

11.65 

4.31 

3.53 

0.21 

1.7 

3.1 

0.4 

0.51 

7.0 

0.20 

0.62 

1.68 

12.13 

8.95 

11.73 

0.26 

8.6 

6.3 

3.4 

0.47 

7.0 

0.30 

0.62 

1.52 

12.75 

15.34 

18.98 

0.28 

17.3 

8.5 

8.8 

0.40 

7.0 

0.40 

0.63 

2.12 

18.27 

0.94 

5.10 

0.28 

4.6 

4.4 

2.3 

0.45 

7.0 

0.50 

0.63 

1.65 

15.39 

2.06 

16.80 

0.33 

20.4 

8.6 

16.1 

0.50 

7.0 

0.75 

0.64 

0.70 

6.95 

12.43 

30.98 

0.32 

36.6 

11.4 

25.3 

0.37 


TABLE 1 

Summary of simulations followed through merger. 


For VpjM > 7.5 only the first fly-by encounter was modeled, hence no information related to disrupted material is available. The energy and angular 
momentum emitted in GWs for these cases drops with increasing Vp after the first encounter, as expected. To within the estimated 20% truncation 
error inferred from the rp = 6 case resolution study, we see no variation with spin. However, even a small variation in the energy emission at fly-by 
could result in a significant change in the time to the subsequent close encounter in a highly eccentric binary. Thus, higher resolution studies would 
be needed to ascertain the effect of spin on the GW signal for VpjM > 7.5. 

^ADM angular momentum. 

^Total energy emitted in GWs through the r = lOOM surface. 

‘^Total angular momentum emitted in GWs. 

"^Bound rest mass outside the BH ~ 10 ms post-merger in percent of Mq. 

®Unbound rest mass in percent of Mq. 

^Rest-mass averaged asymptotic velocity of unbound material. 

^Kinetic energy of ejecta in units of 10^^ erg. 

^Kilonovae bolometric luminosity in units of 10"^^erg s“^ using Eq. (|^. 

^Specific brightness from ejecta interaction with ISM in units of mJy using Eq. (HI). 

J Remnant BH dimensionless spin. 

^Values in parentheses are Richardson extrapolated values using all three resolutions. 
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on the impact parameter and NS spin, these mergers can 
produce accretion disks of up to a tenth of a solar mass. 

We find that near the critical impact parameter the NS 
spin influences the orbital dynamics to a sufficient extent 
to affect whether a merger or fly-by occurs, with a cor¬ 
responding large effect on the GW emission. At a first 
glance this variability might seem exceedingly rare, re¬ 
quiring a finely tuned impact parameter. However, since 
the primary source of this sensitivity to binary parame¬ 
ters arises because the pericenter gets close to the region 
of unstable orbits, which exists for all eccentricities, not 
merely the initially hyperbolic case considered here, one 
can speculate that the last few encounters for any case 
where non-negligible orbital eccentricity remains will be 
subject to this sensitivity. Likewise, the variability as¬ 
sociated with EM counterparts could also be present for 
a larger range of initial impact parameters. Future sim¬ 
ulations of multi-burst events will be needed to address 
this speculation. At the other end of the spectrum, some 
fraction of dynamical-capture binaries that form at larger 
initial separations will circularize prior to merger due to 
GW emission; the results found here thus also motivate 
the study of quasicircular mergers involving millisecond 
NSs. 

We have shown it is important to include spin to under¬ 
stand the full range of possible EM and GW outcomes in 
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eccentric mergers. However, whether it will be possible 
to perform parameter estimation from a putative multi¬ 
messenger event is a different question. Certainly in a 
single burst event the degeneracies will be too strong to, 
for example, identify NS spin as the sole reason for an 
unusually bright counterpart. Multi-burst events can in 
principle lift much of the degeneracy, as information in 
the timing of the bursts could significantly narrow the 
parameters of the progenitor binary. The range of viable 
NS EOSs, NS spin directions, and BH spins needs to be 
simulated, both to determine how these parameters af¬ 
fect the observable outcomes, and how they add to or 
lift degeneracies. GW detection rates and parameter es¬ 
timation also needs to be investigated within a realistic 
data analysis framework including detector noise. All of 
these problems we leave for future studies. We also plan 
to study the effect of spin in dynamical-capture NS-NS 
mergers. 
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